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We present a study of the charmless semileptonic B-meson decays B + — > lu£ + v and B + — > nl + v. 
The analysis is based on 383 million BB pairs recorded at the T(45) resonance with the BABAR 
detector. The u> mesons are reconstructed in the channel u — > 7r + 7r~7r° and the n mesons in the 
channels r\ — » 7r + 7r~7r and r\ — > 77. We measure the branching fractions B(B + — ► ujI + v) = 
(1.14 ± 0.16 s t„t ± 0.08 SHSt ) x 10~ 4 and B{B+ -> V £+v) = (0.31 ± 0.06 stat ± 0.08 svst ) x 10~ 4 . 
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Measurements of branching fractions of charmless 
scmilcptonic B decays can be used to determine the mag- 
nitude of the Cabibbo-Kobayashi-Maskawa matrix [l[ el- 
ement V u b and thus provide an important constraint on 
the Unitarity Triangle. Studies of exclusive decays allow 
for more stringent kinematic constraints and better back- 
ground suppression than inclusive measurements. How- 
ever, the predictions for exclusive decay rates depend on 
calculations of hadronic form factors and are thus af- 
fected by theoretical uncertainties different from those 
involved in inclusive decays. The description of semilcp- 
tonic decays requires one or three form factors for final 
states with a pseudoscalar or a vector meson, respec- 
tively, if lepton masses are neglected. Currently, the most 
precise determination of \V u b\ with exclusive decays, both 
experimentally and theoretically, comes from a measure- 
ment of B — ► -kIv decays 0. It is important to study 
other semileptonic final states with a pseudoscalar or a 
vector meson to perform further tests of theoretical calcu- 
lations and to improve the knowledge of the composition 
of charmless semileptonic decays. 

In this paper, we present measurements of the branch- 
ing fractions B(B + — * ujt + v) and B(B + — > r)£ + v), where 
I = e, \i and charge-conjugate modes are included implic- 
itly These decays have previously been studied by the 
CLEO and BaBAR [IS] collaborations (B+ -> r)l+v) 
and by the Belle @] collaboration (B+ -> loI + v). The 
w meson is reconstructed in its decay to three pions 
(B(u -* tt+tt-tt ) = (89.2 ± 0.7)% 0), while for the 
77 meson the decays to three pions and to two photons 
(B(r) -> tt+tt-tt ) = (22.68 ± 0.35)%, B(rj -> 77) = 
(39.39 ± 0.24)% @) are used. In contrast to earlier 
B + — > i]i + v analyses from BABAR 0, Hj], the second B 
meson in the event is not reconstructed; this yields a 
much larger candidate sample. 

The results presented here are based on a data sample 
of 383 million BB pairs recorded with the BABAR detec- 
tor Q at the PEP-II asymmetric-energy e + e~ storage 
rings at the Stanford Linear Accelerator Center (SLAC). 
The data correspond to an integrated luminosity of 
347 fb _1 collected at the T(45) resonance. In addition, 
35 ib^ 1 of data collected about 40 McV below the res- 
onance (off-resonance) are used for background studies. 
Simulated BB events arc used to estimate signal efficien- 
cies and shapes of signal and background distributions. 
Charmless semileptonic decays are simulated as a mix- 
ture of three-body decays B — > X u iv {X u = tt, 77, 77', p, w) 
and have been reweighted according to the latest form- 
factor calculations from light-cone sum rules [t| [Hi [ll[ . 
Decays to non-resonant hadronic states X u with masses 
mx u > 2-171^ are simulated using the differential decay 
rate given in Ref. [12], which produces a smooth mx u 
spectrum. The GEANT4 package [l3| is used to model the 
BABAR detector response. 

The reconstruction of the signal decays B + — > oj(. + v 
and B + — > i]£ + v requires the identification of a charged 



lepton (e or /i) and the reconstruction of an u> or 
77 meson. The two dominant sources of background 
are semileptonic decays with a charm meson in the fi- 
nal state, B -> X c iv (X c = D,D*,D**,D^tv), and 
e + e~ — > qq (q — u, d, s, c) continuum events. Other 
backgrounds include charmless semileptonic decays that 
are not analyzed as signal and BB events with lepton 
candidates from secondary decays or from misidentifica- 
tion of hadrons as leptons. The center-of-mass momen- 
tum of the lepton is restricted to > 1.6 (1.0) GeV 
[14] ] for the oj (77) final state. This requirement sig- 
nificantly reduces those background events that have 
hadrons misidentified as leptons and rejects a large frac- 
tion of leptons from secondary decays or photon con- 
versions. For the reconstruction of the to or 77 meson, 
charged (neutral) pions are required to have a momen- 
tum in the laboratory frame above 200 (400) MeV to 
reduce combinatorial background. Neutral pion candi- 
dates are formed from two photons with energies above 
100 MeV and an invariant mass, m 77 , in the range 
100 < m 77 < 160 MeV. A three-pion system is accepted 
as an lu (77) candidate if its invariant mass, m^, is in 
the range 760 < m 3w < 806 MeV for w candidates and 
540 < 777,3^ < 555 McV for 77 candidates. The 77 meson is 
also reconstructed via its decay into two photons, each 
with an energy above 50 MeV, with a two-photon invari- 
ant mass in the range 520 < m 77 < 570 MeV. To reduce 
the combinatorial background, two-photon combinations 
are rejected as possible 77 candidates if either of the pho- 
tons can be combined with any other photon of the event 
to form a system with an invariant mass close to the tt° 
mass, 110 < m 77 < 160 MeV. 

Event-shape variables that are sensitive to the topo- 
logical differences between jet-like continuum events and 
more spherical BB events are used to suppress back- 
grounds from e + e~ —> qq and other QED processes. The 
normalized second Fox- Wolfram moment R2 [H| is re- 
quired to be less than 0.5 and a loose requirement on the 
second Legendre moment Li of L2 < 3.0 GeV is im- 
posed. In addition, the event must contain at least four 
charged tracks. 

The charged lepton is combined with an ui{rj) candidate 
to form a so-called Y pseudo-particle candidate, whose 
four-momentum is defined as the sum of the correspond- 
ing lepton and hadron four-momenta. All charged tracks 
belonging to the Y are fit to a common vertex. This ver- 
tex fit must yield a x 2 probability of at least 0.1%. Multi- 
ple Y candidates per event are possible and all candidates 
are retained. The Y multiplicity is well described by the 
Monte-Carlo simulation. About 96% (98%) of simulated 
B + — > u)t + v (B + — > rji + v) signal events and more than 
90% of all selected data events, which include both signal 
and background, contain only one Y candidate. 

The momentum of the candidate neutrino is calculated 
from the difference between the momenta of the colliding- 
beam particles and the vector sum of the momenta of all 
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detected particles in the event. The energy of the candi- 
date neutrino is obtained as the magnitude of its momen- 
tum, since this is less susceptible to bias from lost parti- 
cles or additional tracks than the missing energy, -E m i ss , of 
the event. The magnitude of the missing-momentum vec- 
tor must be at least 500 MeV. The effect of losses due to 
detector acceptance on the reconstruction of the neutrino 
candidate is reduced by requiring the missing-momentum 
vector in the laboratory frame to point into the polar- 
angle range 0.3 < m i ss < 2.2 rad. If the missing energy 
and momentum in the event come from a single unde- 
tected neutrino and the rest of the event is correctly re- 
constructed, the missing mass, TO m i ss , measured from the 
whole event should be compatible with zero. Because the 
missing-mass resolution varies linearly with the missing 
energy, only events with \m^ iss / (2E miBS )\ < 2.5 GeV are 
selected. 

If the Y candidate originates from a signal decay that 
has been correctly reconstructed, the cosine of the angle 
between the B meson and the Y candidate can be calcu- 
lated as cos6 BY = (2E* B E Y - m| - m Y ) / {2\p B \\p Y \). 
Here m b , E B , p B , my , E Y , p Y refer to the masses, en- 
ergies, and momenta of the B meson and the Y can- 
didate, respectively. In the calculation of cos 9 by, the 
.B-meson energy E B and momentum p B are not mea- 
sured event by event. Instead, E* B = yfs/2 is given by 
the center-of-mass energy of the colliding beam particles, 
y/s, and the magnitude of the B momentum is calculated 
as \p B | = \J E* B — m' B . Signal candidates are required to 
satisfy — 1.2 < cos6>by < 1.1. This requirement was kept 
loose to account for the limited detector resolution and 
photon energy losses. 

To reduce backgrounds without significant loss of sig- 
nal, the momenta of the lepton, p* t , and of the hadron, 
P£j „, that make up a Y candidate are restricted. For 
B + — > ui£ + i/, the momenta are required to satisfy \p^\ > 
1.3 GeV or > 2.0 GeV or \pZ\ + \p* £ \ > 2.65 GeV. In 
the case of B + — > rjl + i/, the conditions \p^\ > 1.3 GeV 
or \p}\ > 2.1 GeV or \p* n \ + \p* t \ > 2.8 GeV have to be 
fulfilled. 

The kinematic consistency of the reconstructed Yv 
system with a signal B decay is verified using 
the two variables AE = (Pb ■ -Pbeam — s /2)/V* 
and toes 



TABLE I: Signal efficiencies, ^signal, and signal-to- 
background ratios, S/B, after the neural- network selection. 



VP 



-Pbeam — (-^bcam j Ph 



where 



PB ■ Pbeam) I "beam 

) is the four-momentum of the 
colliding beam particles and Pb = (Eb,Pb) is the B- 
meson four-momentum computed as the sum of the four- 
momenta of the Y and the v candidates. These vari- 
ables are later used to extract the signal yields in a fit 
to the two-dimensional AE vs. toes distribution. Only 
candidates that fulfill the loose requirement \AE\ < 
0.95 GeV and toes > 5.095 GeV (fit region) are retained. 

At this stage of the selection, the signal-to-background 
ratio, S/B, where S and B denote the expected signal 
and background yields, respectively, is small. It amounts 
to 1.5% for B+ -> ul+v, and 1.8% (1.0%) for B+ -> 









Fit region 


Signal region 










S/B 


S/B 


B+ 


-y lj£ + i> 




1.00 


0.15 


0.46 


B+ 


— > r)l + V, 77 — » 7T + 7T~ 


V 


1.62 


0.10 


0.35 


B+ 


— ► rjl + v, rj — > 77 




4.90 


0.04 


0.15 



T]i + v with rj — ► 7r + 7r~7r°(77). The signal efficiencies for 
the sum of decays with electrons and muons, estimated 
from simulation, are 2.8% for B + — > oj£ + v and 4.0% 
(9.4%) for B+ -> r}t + v with r/ -> tt+ tt~ tt° (77) . 

For further discrimination between the signal and the 
background, a multivariate selection based on neural net- 
works [17[ is used. For each of the three signal channels 
under study, neural networks with two hidden layers (four 
and two neurons, respectively) are applied consecutively 
to separate the signal from the two main backgrounds. 
A first neural network discriminates the signal against 
qq continuum events; a second network is used to further 
distinguish the signal from the B — > X c £v background. 
The neural-network decision is based on the following in- 
put variables: m 2 mlss / (2E mlss ), 6> mlss , cos 9 by, R2, L 2 , 
cos AOthrust , the cosine of the polar-angle difference be- 
tween the thrust axes of the Y candidate and of the rest 
of the event, and cos 9we , the cosine of the lepton "he- 
licity angle" measured in the rest frame of the virtual 
W (calculated using the lepton and neutrino candidates) 
relative to the W direction in the B rest frame. For the 
three-pion final states, the Dalitz amplitude, the magni- 
tude of the vector product of the ir + momentum and the 
7r~ momentum in the cj/77 rest frame, normalized to its 
maximum value, serves as an additional input variable 
to separate u> 7r + 7r~7r° and 77 — > 7r + 7r~7r° decays from 
combinatorial background. 

The training of the neural networks is done using the 
corresponding simulated signal and background samples 
for each of the three signal channels separately. Inde- 
pendent simulated event samples arc used to validate the 
training. Based on Monte-Carlo simulation, the selection 
criterion for each of the output discriminants is chosen 
to maximize the quantity S/ y/S + B. The signal effi- 
ciencies and the S/B ratios after the neural-network se- 
lection are given in Table Q] for the fit region and for the 
signal region delimited by —0.2 < AE < 0.4 GeV and 
m ES > 5.255 GeV. 

For the determination of the signal branching fractions, 
the AE vs. toes distributions of the simulated signal 
and backgrounds are fit to the data distribution for the 
three signal channels independently. The fits are based 
on an extended binned maximum-likelihood method (l8j 
and take statistical fluctuations of both the data and the 
Monte-Carlo samples into account. The binning of the 
AE vs. Toes distributions used in the fits contains a total 
of 50 bins with smaller sizes in the signal region to re- 
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solve the signal shape and larger sizes in the part of the 
fit region outside of the signal region to determine the 
background normalizations from data. The shapes of the 
signal and background distributions are taken from sim- 
ulation. The fits determine the relative fractions of the 
signal and some of the background samples in the data. 

The free parameters of the fits are the normalizations 
of the signal and the B — > XJv background, and for the 
B + — > cu£ + v channel also the overall normalization of the 
continuum background. The B — ► XJv normalization 
is left free to account for a slight discrepancy between 
the B — > XJv yields in data and Monte-Carlo simu- 
lation. The relative contributions of events with elec- 
trons or muons to the continuum background have been 
determined using off-resonance data. Compared to the 
B + — > uji + v channel, the B + — > r]l + v channels suffer 
from a larger continuum background and the fit shows 
larger correlations between signal and this background 
component. The normalization of the continuum back- 
ground in the B + — ► r]£ + v channel is therefore taken 
from off-resonance data and is not varied in the fit. All 
other background distributions, mainly other B — > X u £v 
decays and B decays with secondary or misidentified lep- 
tons, are fixed to their Monte-Carlo predictions. The fit 
procedure has been validated with simulated signal and 
background samples. 

The resulting signal yields and branching fractions 
for the three signal channels are presented in Table [Til 
The scale factors between the B — ► X c lv background 
yields predicted by the simulation and the values deter- 
mined by the fits are 1.06 ± 0.07 for B + — > io£ + v and 
0.96±0.07 (1.12±0.03) for B+ -> rj£+v with r\ -> tt+tt-tt 
(jj — > 77). The correlations between the signal and the 
B — > XJv parameters determined by the fits are 0.08 
for B+ -> ujI+v and -0.60 (-0.46) for B+ -> <ql + v 
with rj — > 7T + 7r _ 7r° (77 — > 77). The correlation be- 
tween the signal and the continuum parameters for the 
B + — > cu£ + v channel is —0.55 and the continuum back- 
ground normalization is adjusted by a factor of 0.89±0.12 
with respect to the normalization obtained from the off- 
resonance data sample. The goodness of fit is evaluated 
using a x 2 -based comparison of the fitted AE vs. jjjes 
distributions of the simulated and data samples and is 
shown in Table [TTJ In addition, the combined branch- 
ing fraction for the two B + — > rj£ + v channels has been 
obtained from a fit to the sum of the AE vs. tties distri- 
butions for r) — > 7T + 7r~7r° and rj — > 77. 

Figure [T] shows the projections of the fitted distribu- 
tions on AE and toes for the three signal channels and 
the combined B + —> rj£ + v channel. For a better vis- 
ibility of the signal, the AE projections are shown for 
mES > 5.255 GeV and the ijies projections are shown 
for -0.2 < AE < 0.4 GeV. 

The systematic errors on the measured branching frac- 
tions arc listed in Table IIIII They arc estimated by 
varying the detection efficiencies or the parameters that 



TABLE II: Signal yields and corresponding branching frac- 
tions as determined by the fits for the three signal channels 
and the combined B + — > rj£ + v channel. The last row shows 
the x 2 P er degree of freedom. 

\B + -> luI + u\ B + -> 7]t + v 

77 — ► 7r + 7r - 7r° r) —* 77 combined 
N signal 802 ± 113 127 ± 42 459 ± 98 554 ± 105 
B(W~ 5 ) 11.4 ± 1.6 4.36 ± 1.43 3.01 ± 0.64 3.11 ± 0.59 
X 2 /d.o.f. 36.0/47 59.9/48 43.2/48 49.7/48 



impact the modeling of the signal and the background 
processes within their uncertainties. The complete anal- 
ysis is then repeated and the differences in the result- 
ing branching fractions are taken as the systematic er- 
rors. The total systematic error is obtained by adding in 
quadrature all listed contributions. 

Uncertainties due to the reconstruction of charged par- 
ticles and photons are evaluated by varying their recon- 
struction efficiencies and the energy depositions of pho- 
tons in the simulation. The neutrino reconstruction is 
affected by background with long-lived K®, which often 
escape detection and contribute to the measured missing 
momentum of the event. The uncertainty arising from 
K® L production and interactions is estimated by varying 
their production rate as well as their detection efficiency 
and energy deposition in the simulation. For lepton iden- 
tification, relative uncertainties of 1.4% and 3% are used 
for electrons and muons, respectively. A 3% uncertainty 
is assigned to the ir° /rj — > 77 reconstruction efficiency. 

The uncertainty due to the B — > XJv background is 
evaluated by varying the B — ► D/D* /D**£v branching 
fractions Q and the B — > D* form factors Prior to 
the neural-network selection, the background level is high 
and discrepancies between data and Monte-Carlo distri- 
butions are observed at roughly the 10% level for some of 
the neural-network input variables. To estimate the effect 
of these discrepancies on the measured branching frac- 
tions, the dominant background component (B — > XJLv) 
is reweighted to match the data. The weights are de- 
termined from a B — > X c £^-enhanced sample which is 
obtained by selecting only events that are otherwise re- 
jected by the B — ► XJv neural-network selection and 
keeping all other selection criteria unchanged. 

For the B — > X u lv background, the branching frac- 
tions of the exclusive decays that are not analyzed as 
signal are varied within their uncertainties [20(. The 
non-resonant part is varied within the range allowed by 
the uncertainty of the total B — > X u £v branching frac- 
tion [2(J. The uncertainty due to the normalization of 
the continuum background has been determined with 
off-resonance data for events with electrons or muons 
separately. Since the overall normalization of the con- 
tinuum background is adjusted for B + — ► uj£ + v in the 
fit, the resulting error in this channel is smaller than 
for the B + — > i]£ + v channels. For the normalization 
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FIG. 1: (color online) Projected AE distributions for /ties > 5.255 GeV (top), and tues distributions for —0.2 < AE < 0.4 GeV 
(bottom). From left to right: B + — * toi + v channel, B + — > rj£ + v channel with 77 — > -k + tv~-k , B + — » rj£ + v channel with 77 — » 77, 
and combined B + — > rjl + v channel. The error bars represent the statistical uncertainties on the data. The histograms show 
simulated distributions for signal (white), B — > X c £u decays (light shaded/yellow), gg-continuum (dark shaded/blue) and all 
other backgrounds (hatched) and have been summed up. The distributions of the simulated signal and B — > X c iv background 
(and the qq background for B + — > uj£ + v) have been scaled to the results of the fits. 



of secondary-lepton background, an uncertainty of 6-8%, 
depending on the signal channel, has been estimated from 
a detailed study of the composition of this background. 

Uncertainties in the modeling of signal decays due to 
the imperfect knowledge of the form factors affect the 
shapes of kinematic spectra and thus the acceptances of 
signal decays. The errors on the measured branching 
fractions are estimated by varying the parameters of the 
form-factor calculations withing their uncertainties 

a 

The branching fractions of the decays uj/rj — ► tt + tt tt° 
and 77 — ► 77 are also varied within their uncertain- 
ties [7J . The uncertainty on the number of produced B 
mesons is 1.1% [2l|. The uncertainty on the ratio of 
the T(45) -> B+B~ and T{AS)_^ B°W branching frac- 
tions, /-| //oo = 1.065±0.026 [20||, is taken into account. 

The total systematic errors on the measured branching 
fractions are 7.2% and 25.1% for the B + — > ujI + v and the 
combined B + — ► rj£ + v channels, respectively. 

In summary, we have measured the branching fractions 
of B + — > loI + v and B + — > rj£ + v decays to be 

B{B + ^ujI + v) = (1.14 ± 0.16 ± 0.08) x 10" 4 , 
B{B+^i^ + v) = (0.31 ± 0.06 ± 0.08) x 10" 4 , 

where the errors are statistical (data and simulation) and 
systematic, respectively. 

The B + — > rjl + v and B + — > lo£ + v measurements pre- 
sented here significantly improve the current knowledge 



of these decays. The B + — > r]£ + v result is compati- 
ble with an earlier measurement by BABAR 0| based on 
events tagged by a semileptonic decay of the second B 
meson, B(B+ -» r]t+u) = (0.64 ± 0.20 stat ± 0m syst ) x 
10 -4 . The two analyses are statistically independent 
and complement each other. The analysis presented 
here is statistically more precise but has larger system- 
atic uncertainties, as expected for an untagged measure- 
ment. We combine the two BABAR results and obtain 
B(B+ -> r)t+v) = (0.37±0.06 stat ±0.07 syst ) x lO" 4 . 
The B + — > uot + v branching-fraction measurement is the 
first with a significance of more than five standard devi- 
ations. It represents an improvement by a factor of three 
over the only earlier measurement by Belle @. The im- 
proved measurements of B + — * uj£ + iy and B + — * r]t + v 
decays are important ingredients to the determination 
of the composition of the inclusive charmless semilep- 
tonic decay rate. The size of the data samples is not yet 
sufficient to perform a measurement in intervals of the 
momentum transfer, q 2 , of the decay, which would be 
necessary to determine \V u b\ with an adequate precision. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BABAR. The collaborat- 
ing institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
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TABLE III: Relative systematic errors of the branching fractions B(B + — > lu£ + u) and B{B + — > rjl + v). For the B + — > ?7^ + f 
channel, the systematic errors for the three-pion and the two-photon final states as well as for the combined result are shown. 
The total error in each column is the sum in quadrature of all listed contributions. 



Error source 


8B{B + tul + v) (%) 


5B{B + 


-> rjl+v) 


(%) 






r\ — > TV + TT~ir° 


77 -» 77 


combined 


Tracking efficiency 


1.9 


4.9 


4.2 


4.6 


Photon reconstruction 


2.1 


1.8 


9.1 


8.6 


Kl production and interactions 


2.6 


4.8 


3.1 


1.9 


Lepton identification 


1.9 


3.3 


6.9 


6.3 


Tr°/rj identification 


3.8 


6.9 


13.3 


12.2 


Neural-net input variables 


0.6 


0.8 


5.9 


6.1 


D* form factors 


0.4 


1.0 


0.9 


1.0 


B(B -> XJu) 


2.1 


5.5 


7.6 


8.0 


B(B -> X a tv) 


2.8 


4.4 


9.8 


8.6 


Secondary leptons 


0.3 


0.3 


0.6 


0.5 


Continuum scaling 


0.7 


15.8 


10.4 


12.7 


Signal form factor(s) 


1.8 


5.9 


0.3 


1.3 


B(u>/r] — > 7r + 7r _ 7r°), B(r\ — > 77) 


0.8 


1.5 


0.6 


1.2 


N B b 


1.1 


1.1 


1.1 


1.1 


/+-//00 


1.2 


1.4 


0.9 


1.0 


Total systematic error 


7.2 


21.1 


25.2 


25.1 
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